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Abstract. We establish rigorously the existence of a three-parameter family of self-similar, 
globally bounded, and continuous weak solutions in two space dimensions to the compress- 
ible Euler equations with axisymmetry for 7-law poly tropic gases with 1 < 7 < 2. The 
initial data of these solutions have constant densities and outward-swirling velocities. We 
use the axisymmetry and self-similarity assumptions to reduce the equations to a system 
of three ordinary differential equations, from which we obtain detailed structures of solu- 
tions besides their existence. These solutions exhibit familiar structures seen in hurricanes 
and tornadoes. They all have finite local energy and vorticity with well-defined initial and 
boundary values. 



1. Introduction. 

We are interested in finding some solutions to the initial value problem for the two- 
dimensional compressible Euler equations. Our approach is to try to generalize to two 
dimensions some of the results on Riemann problems for the one-dimensional compressible 
Euler equations. One natural generalization is to consider initial data which consist of 
four constant states, or any finite number of constant states [13-15]. The well-known 
configurations of regular and Mach reflections are special cases of such a generalization. 
However, no rigorous proof of existence of any nontrivial solutions to these problems has 
been established. 

A more complete generalization is to consider initial data which depend only on the 
polar angle in the two-dimensional space of positions. This generalization looks certainly 
more formidable, if one hopes to solve it all. However, it now contains a special three- 
parameter family of data, namely, the axisymmetric initial data, which allows us to assume 
axisymmetry of the solutions and reduce the initial value problem of the partial differential 
equations to a boundary value problem of a system of ordinary differential equations. 
With this symmetry we were able to construct in [16] a two-parameter family of selfsimilar 
solutions corresponding to pure rotational initial data. In paper [17], we gave rigorous 
proofs for these solutions and constructed an additonal one-parameter family of solutions 
for the square polytropic cases. This additional parameter allows the initial flows to swirl 
outward. In the current paper we take on the sub-square cases: we construct rigorously a 
three-parameter family of solutions. 

Our main task here is to study the resulting boundary value problem of the non- 
autonomous system of three ordinary differential equations in the cases 1 < 7 < 2. Singu- 
larity (stationary) points of the system consist of two-dimensional manifolds in the four- 
dimensional phase space. These singularity points correspond to surfaces of characteristics, 
where solutions may not be differentiable in the physical space and time. Linearizations at 
some stationary points yield no structure for the solutions due to the high order degeneracy 
of the system at these points. The existence of connecting orbits follows from high-order 
local analysis at stationary points and global analysis of various invariant regions of the 
system. A typical global solution may consist of as many as three nontrivial connecting 
orbits chained together continuously. The complete construction of the three-parameter 
family of solutions is described in the conclusions section, section 8, at the end of the 
paper. 
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We find that our solutions capture qualitatively some typical properties of swirling 
flows such as hurricanes and tornadoes. As an illustration, we plot a typical solution in 
Fig. 1.1-2. These figures are placed out of the main flow of the text due to their large 
sizes. See Fig. 8.1-2 for more precise description on the eye and wall regions. We refer the 
interested reader to paper [16] for explicit presentations of spiralling particle trajectories. 

We point out that our initial data are restricted to the set of data with nonnegative 
radial velocities (swirling outward), in addition to axisymmetry and radial symmetry. 
Shock waves may arise if the initial radial velocities are allowed to be negative (swirling 
inward). Due to complications of the global dynamics of the ordinary differential equations, 
we have only studied the special case 7 = 2 in a companion paper [17]. In this paper, we 
study the case 1 < 7 < 2. The case 7 > 2 will be considered in a forthcoming paper. For 
some partial results on the case 7 > 2, we refer the reader to our paper [16]. Preliminary 
steps (Sect. 1-5) are reproduced here from [17] for readers' convenience. 

There is a great deal of related work, from which we mention only the most closely 
related. For general existence of weak solutions with axisymmetry for the 2-D compressible 
Euler equations outside a core region, we refer the reader to the recent work of Chen and 
Glimm [5]. For some explicit solutions of the compressible Euler equations with spherical 
symmetry but without swirls, see Courant and Friedrichs [4]. For viscous swirling motions, 
we refer the reader to Bellamy-Knights [1], Colonius, Lele, and Moin [3], Mack [9], Powell 
[11], and Serrin [12]. 



2. The problem. 

We consider the two-dimensional compressible and polytropic Euler equations 

Pt + (pu) x + (pv) y = 0, 

(pu) t + (pu 2 + p) x + (puv)y = 0, (2.1) 

(pv) t + (puv) x + (pv 2 +p) y = 0, 

where p = p(p) is a given increasing function of p. Global existence of weak solutions to 
its initial value problem is open. Attempts have been made through considering special 
situations such as the diffraction of a planar shock at a wedge or generalized Riemann 
problems with four different initial constant states. 

Here we consider a situation which involves swirling motions. We impose axisymmetry 
to the system. That is, we assume that our solutions (u, v, p) have the property 

p(t,r,6) = p(t,r,0) 

fu(t,r,9)\ fcosO -sm6\ f u(t,r,0)\ , , 

\v(t,r,e) J \sin9 cos 9 J \v(t,r,0) J { ' 

for all t > 0, 9 G K. and r > 0, where (r, 9) is the polar coordinates of the (x, y) plane. 
With this symmetry, system (2.1) can be reduced for continuous solutions to 

/ N P u 

p t + {pu)r H =0, 
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2 

u t + uu r + — — — = 0, (2.3) 
p r 

uv 

v t + uv r H = 0, 

r 

where p = p(£, r, 0), etc. Note now that u and v in (2.3) represent the radial and pure 
rotational (a.k.a. tangential ) velocities in the flow, respectively. 

We limit ourselves to a Riemann-type of initial data; that is, we require the initial 
data to be independent of the radial variable r > 0: 

(p(0, r, 9), u(0, r, 0),v(0, r, 9)) = (p (9), u (9), v (9)). (2.4) 

We remark in passing that the initial value problem of system (2.1) with the type of data 
in (2.4) seems to be more in the spirit to be called the two-dimensional Riemann problem 
for (2.1). It degenerates to the classical Riemann problem for the one-dimensional case, it 
is simple, and yet general enough to contain important waves such as swirling motions as 
well as shock and rarefaction waves and slip lines (surfaces). 

When the axisymmetry condition (2.2) is imposed onto (2.4), we find that our data 
are limited to be 

u(0,r,9) = uq cos 9 — vo sin 9, 

v (0,r,6>) = uo sin6> + vq cos6>, (2.5) 
p(0,r,9) = po, 

where p > 0, (u , v ) G M 2 are arbitrary constants. Hence our data for system (2.3) are 

p{0,r,0) = p o , u{0,r,0) = u o , v{0,r,0) = v o . (2.6) 

Since the problem (2.3) (2.6) is invariant under self-similar transformations, we look 
for self-similar solutions (p, it, v) which depend only on £ = r/t. We thus have the following 
boundary value problem of a system of ordinary differential equations: 

pQ IE uv 

Pr=—r, ^ = _ X' v r=— { 2 ' 7 

r A r A ryr — u) 



where 





lim (p,u,v) = 


(Po, uo- 




A = 


c 2 - (u - r) 2 , 






e = 


v 2 — u(r — u), 






E = 


(r — u)Q — uA = 


v 2 (r — 


u)- 


c = 









(2.8) 



•uc 2 , 
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and r is used in place of £ for notational convenience. We note that c is the sound speed. 

We will construct global continuous solutions or establish their existence for problem 
(2.7-8) for any p > 0, uq > 0, vq € R and p(p) = A 2 /0 7 where 1 < 7 < 2 and A 2 > is a 
constant. We shall be mainly concerned with the case 1 < 7 < 2 in Sections 4-7, but we 
will mention the differences for 7 = 1 in Section 8. We will also assume vq > because 
the case v < can be transformed by v — > — 1> to the case vq > 0. It can be verified that 
these solutions are also global continuous solutions to the original Euler equations. 



3. Far-field solutions. 

We show that problem (2.7-8) has a local solution near r = +00 for any datum 
(po,uo,vo) with po > 0. Let s = K Then (2.7-8) can be written as 



dp 


p[u(l — us) 


— 


ds 


S 2 C 2 - (1 - 


- us) 2 


du 


sue 2 — v 2 (l 


— us) 


ds 


S 2 C 2 - (1 - 


■ us) 2 


dv 


uv 




ds 


1 — us 




(p, u,v 


)\a=0 = (PO, 


U , V 



(3.1) 



Problem (3.1-2) is a classically well-posed problem which has a unique local solution for 
any initial datum with p > 0. 
We find that v = and 

sv 2 - u(l - us) = (3.3) 
are invariant surfaces in the four-dimensional (p,u,v,s) space. It can be verified that 

d r 2 m (1 _ 3SW)SC 2 + u(l - SU) 2 2 

— [SV -U(1-US)\ = T — — -^-\SV -U[\-US)\. 

ds 1 — su )\s z c z — 1 — su Y\ 



Explicit solutions. In the special case uq = 0, we find in [16] from the invariant 
surface (3.3) a set of explicit solutions near r = +00 : 



p = p , u=f, v = ^^r 2 -vl r>r* (3.4) 
r* = \ (V c o + 4«g + c o) , (3.5) 



where 



and 

co = Vp'(Po)- 
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We comment that the functions in (3.4) are actually defined for all r > vo, but we cannot 
use them up to vq with absolute certainty since (3.1) has a singularity at the point r = r* 
on the curve (3.4). In fact, we find that the position r, the radial velocity u, and the sound 
speed c at r* along (3.4) have the relation 

r* = u(r*) + c , 

which is to say that r* is the radial characteristic speed; that is, distance that a small 
disturbance generated from the origin at time zero can travel radially in time t = 1. Thus, 
r* in (3.5) is a characteristic position. 

We point out that the Mach number along any solution in (3.4) is a constant: 

M = \/ u 2 + v 2 /c = vq/c . 

For those who are interested in the pseudo Mach number defined in the self-similar coor- 
dinates (£, rf) as 

M s EE y/( U -0 2 + (v-V) 2 /c, 

we find that 

M s = ^ r 2 -v 2 /c 
along a solution of (3.4). At the end point r = r* , we find 

M s (r*)= ^M 2 + l/4 + l/4^ ' >1, 

where 

M ee v /c 

is the Mach number for the initial data in this case. Pseudo Mach number determines 
the type (hyperbolic or other) of the 2-D Euler equations in the self-similar plane (£, rf): 
hyperbolic if M s > 1, other types if M < 1, see [15] for example. So the solutions (3.4) 
are pseudo supersonic and in hyperbolic regions. 



4. Intermediate field equations. 

We can simplify system (3.1) by assuming the polytropic relation 

p(p) = A 2 p^ (4.1) 
for some A 2 > and 7 > 1, and introducing the variables 

I = su, J = sv, K = sc. (4-2) 
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Then system (3.1) can be written into the form 

dl 2IK 2 -(1-I)[J 2 + 1(1-1)] 



s 



ds K 2 -(l-I) 2 



dJ T 1-2J 

s ds~ = J T^T < 4 - 3 > 

dK_ _ K 2K 2 - 2(1 - I) 2 - (7 - 1) [J 2 - 1(1 - I)] 
S ~ds~ ~ ~2 K 2 - (l-I) 2 ' 

Corresponding to the initial data (3.2), we shall look for solutions of (4.3) with the following 
initial condition: 

(/, J,K) ~ s(u ,v ,co) (4.4) 

as s ^ 0+. We note that (4.3) is now autonomous for 7, J, K with respect to the new 
variable s' = In s. 

The invariant surfaces of (4.3) are the surface J = 0, the surface K = 0, and the 
surface 

H = J 2 - 1(1 - I) = (4.5) 

since 

d (1 - 2I)[2K 2 - 3(1 - I) 2 ] 
ds (1 - I)[K 2 - (1 - I) 2 ] ' 

The explicit solutions (3.4) are all in the invariant surface H = 0. We postpone the 
description of these solutions in the variables (I, J, K) to Sect. 6 when we have a more 
complete picture of all solutions. 

Scaling symmetry. System (4.3) is invariant under the coordinate transformation 
s — > as for any constant a > 0. In particular, we can take a = 1/cq. Thus we may assume 
that po > is such that Co = \Jp'(po) = 1. Hence the structure of any solution of problem 
(4.3-4) will depend only on the dimensionless quantities uq/cq and vq/cq. 



5. Solutions without swirls. 



Let us first determine the distribution of integral curves on the invariant surface J = 0. 
Assume v = 0. We look for solutions to problem (4.3-4) with J = 0. Hence we have 
a subsystem for (J, K): 

(5.1) 



di _ T 2K*-{i-iy 

b ds ~ 1 K 2 -(l-I) 2 



,dK 



K K 2 -(1-If + ^ I(1-I) 



K 2 -{i-iy 

Introducing a new parameter r, we can rewrite (5.1-2) as 



' dI =I[(1-I) 2 -2K 2 } 



dT 



^ = K[(l- I) 2 - K 2 - ^1(1- I)] 
[f T =s[(l-I) 2 -K 2 ]. 
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(5.2) 

(5.3) 
(5.4) 
(5.5) 



Note that (5.3-4) form an autonomous subsystem. If uq also vanishes, then we have a 
trivial solution p = po? u = v = 0. 

Suppose uo > 0. It can be verified that our far-field solutions starting at s = 0+ enter 
the region Q C M 2 in the (I, K) phase space given by 

{/ > 0, K > 0, 
a = (1 _ j)2 _ ^2 _ 3-1 j(l _ j) > in 0</<± 
b = (1 - I) 2 - 2K 2 > in i < / < 1 

See Figure 5.1. 




Fig. 5. 1. Phase portrait of solutions without swirls. 

We show that solutions starting in the closure Q will not leave O as s increases. 
Notice first that s > is an increasing function of r in O by equation (5.5) so we will show 
that solutions of the two equations (5.3-4) starting in the closure O will not leave Q as r 
increases. The stationary points of (5.3-4) in are the points (I, K) = (0, 0), (0, 1), (1, 0), 

and 75 (l _ ^)) • Tne axis if = 0in0</<l, and the axis / = in < K < 1 are 
trivial solutions. On the boundary 6 = 0, i.e., 

K=-L(l-J), I<I<1, (5.6) 

we find that j£ = 0, ^ < 0. So solutions enter O on (5.6). Finally on the boundary 
a = 0, i.e., 

K 2 = (l-/) 2 - l^ij(l-J), 0</<-, (5.7) 
2 7 
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we find that ^ < 0, and ^ = 0. So solutions enter O on (5.7) also. 

We conclude without showing further details that some solutions in O will go to 
the point (1,0), while others go to the point (0,1), with exactly one integral curve (the 

heteroclinic orbit) going to the point (^"^ffl - 7 ) )' 




Fig. 5.2. 



5.1. The transitional solutions. To determine what data {uq, Po) yield the transitional 
solution leading to the stationary point ^1 — we first eliminate r from (5.3-4) 

by division and introduce L = (K/I) 2 to find 

dI-~ L (1-IY-2PL ' 0</<1 ( 5 - 8 ) 
L = Mq 2 at 1 = (5.9) 

where 

M = u /c 

denotes the Mach number of the initial states (-u ,0, p )- Problem (5.8-9) is well-posed 
for any M > 0. The transitional solution goes from the point (I, L) = (0, M ~ 2 ) to 

(7' 1(7 — I) 2 )- However, there does not seem to have an explicit formula for the value 
Mq which yield the transitional solution. We give an estimate instead. For convenience 
we let Mh('y) denote the value of the initial Mach number for the transitional (critical, 
heteroclinic) solution. It can be seen that M = is an upper bound for the transitional 

Mach number Mh("f). In fact, any solution L(I) of (5.8-9) with datum M > starts as 

a decreasing function of / > till / = K In the interval I e ^, lj , the solution remains 
below the two curves on which the numerator and demoninator of the right-hand side 
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of (5.8) vanish respectively, therefore remains decreasing until the final stationary point 
(J, L) = (1, 0), see Figure 5.2, where d = (7 — l) 2 /2, / = M^ 2 (7), and N and D are where 
the numerator and denominator of the right-hand side of (5.8) vanish respectively. 

This transitional solution yields a one-parameter family of smooth solutions in terms 
of (r,u,v,p). We see from equation (5.5) that Ins approaches infinity as the solutions 

approach the point (I,K) = ^1 - ^jj, since (1 - I) 2 - K 2 7^ at the point. So 

s — > +00 and thus r — > 0+. Also the solutions have the asymptotics 

u(r) = — r, c(r) = —= ( 1 ) r 

7 V2V 7/ 

as r — > 0+. These global transitional solutions are similar and one is sketched in Figure 
5.3. 




Fig. 5.3. A zero-swirl transitional solution. 



The Mach number at r = of the transitional solution is M = v / 2/(7 — 1). The 
corresponding pseudo Mach number is M s = 

5.2. Finiteness of the parameter s. We show next that the parameter s approaches 
finite numbers when solutions of (5.3-5) approach either the points (I, K) = (1, 0) or (0, 1). 

We can linearize the two equations (5.3-4) at (I, K) = (0, 1) to find 



iL — _ t 

dr ~ ± 



(5.10) 

diEzll = _7±3 / _ 2(K _ 1) . 

The eigenvlaues of (5.10) are Ai = —1, A 2 = —2. So solutions of (5.3-4) near (0, 1) approach 
(0, 1) exponentially as r — > +00. From equation (5.5), we find 



lni. = f Ul-I) 2 -K 2 ]dr 

S Jro 
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for some constants sq > and To- So s approaches a finite number as r — > oo since 
(1 — J) 2 — i^ 2 approaches zero exponentially. 

Linearization at the point (I, K) = (1, 0) of the two equations (5.3-4) yield the trivial 
system of zero right-hand side. We need a different approach. We show that solutions near 
(1,0) will enter (1,0) in the sector bounded by K — and the line 

K = a(l-I) (5.11) 

for some < a < 4=, such that a 2 + > 0. See Figure 5.4. 



K 




(1,0) 



Fig. 5.4. The parameter s is finite near (1,0). 



In fact, a vector in the normal direction of (5.11) is (a, 1). We calculate the inner 
product of the vector field of (5.3-4) with the direction (a, 1) to find 



^(I,K).(a,l) = -a 



(l-a 2 )(l-J) 



{I -If. 



(5.12) 



The expression (5.12) is negative if I is close to 1 and a is such that a 2 + > 0. So we 
conclude that every solution that ends at (1,0) will be such that 



K < a(l-I) 

near / = 1 for some a < since 7 > 1. 

Now we look at the equation (5.1) and use (5.13) to find 

dl / l-2a 2 
ds 1 — a 2 



(5.13) 



11 



when I is close to 1. Thus 

dl 

s— > a positive constant 
as 

when / is near 1. Therefore s is finite at the end (1,0) since / is bounded and the geometric 
integral f^°° ^ds diverges. 

5.3. The global construction. We are now ready to construct global solutions for (5.1- 
2). For each solution ending at (I,K) = (0, 1), we continue the solution by the constant 
state (u, v, p) = (0, 0, p**) where p** is the value of p at the ending point. These are 
continuous extensions. The relation between the density p** and the ending value s** is 

p '(p**) = s **- 2 . 

Let 

c** = VpIpI, 

then 

c = r 

when r** is the radius of the circle of the constant state in the physical plane t = 1. So 
the core region has a constant density which expands at its sound speed. The Mach and 
pseudo Mach numbers at the edge of the core region are and 1 respectively. In the core, 
the pseudo Mach number M s = r/c** < 1. We therefore have constructed global solutions 
in this case. 

For each solution ending at the point (I, K) = (1, 0), we continue the solution by the 
vacuum state p = 0. We do not need to specify the functions u or v in the vacuum since 
the Euler equations have p as a factor in every term. Each vacuum occupies a circular 
region of radius r** determined by r** = = u** in the physical plane at t = 1, where 
u** is the terminal radial velocity of the fluid at the edge of the vacuum. Asymptotically, 
we find (see also [16]) 

u - r = - 2(7 ~ i 1) (r - r**)(l + o(l)), 

c= (7-l)V^7 (r _ r „ )(1+o(1)) 
7 + 1 

at the edge for 7 < 3. For 7 = 3, the asymptotics is given by 

u-r = -(r -r**) (l + o(l)), 

c=-(r- r**) ln" 1/2 (1 + o(l)). 

And for 7 > 3, it is given by 

u-r = -(r-r**)(l + o(l)), 

c = K (r - r**)"^ (1 + o(l)), K > 0. 
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From these we find that the Mach number M = oo for all 7 > 1 and the pseudo Mach 
number M s = for 7 < 3, and M s = 00 for 7 > 3 at the edge. 

In all, we have constructed global solutions to the reduced system (5.1-2), the special 
case of (4.3-4) with zero swirl. 

6. General solutions in the intermediate field. 

Now consider the case v > and u > for system (4.3). Let fi 3 C M 3 be the set of 
points (I, J, K) satisfying < i" < 1, J > 0, K > 0, 

H = J 2 - 1(1 - I) < 0, 

B=(l- I)[J 2 + 1(1 - I)] - 2IK 2 > if - < I< 1, 

7 

A = 2(1- I) 2 + (7 - 1)[J 2 - 1(1 - I)] - 2K 2 > if < I < -. 

7 

See Figure 6.1. Note that the surface A = intersects the coordinate plane K = on the 
line 

^ = I±l(l-J)(J-^_) 
7-1 7+1 

which lies inside H < 0. It can be verified that all far-field solutions with uq 7^ 0, v$ > 
0, po > enter the region in s > 0. Far- field solutions with uq = enter the side H = 
of 3 , see Sect. 4. We omit these tedious verifications. 

Also let 33 C M 3 be the set of points (I, J, K) satisfying < i" < 1, J>0, < K < 
1 - /, and H < 0. See Fig.6.1. 

We find that it is convenient to introduce a new variable r, as in section 5, to write 
the system (4.3) in the following form 



dl 

dJ 
d7 



(l-I)B (6.1) 



J(l-2I)[(l-iy -K z ] (6.2) 



f = lm-DA (6.3) 

^ = S (l-/)[(l-/) 2 -A' 2 ] (6.4) 

This is an autonomous system for (/, J, K, s), and the first three equations (6.1-3) form an 
autonomous sub-system for (/, J, K ) . 

The stationary points of the system (6.1-3) contained in the closure ^33 are found to 
be given by the edge 

K = 1 — I, J 2 = 1(1 — /), VI G [0, 1] (6.5) 
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and 

(/, J , K )=g,0,-^(l-I)) (6.6) 

in the case 7 7^ 2. For the case 7 = 2, the stationary points of (6.1-3) are given by edge 
(6.5) and the curve 

I = h 1 

2 ' , 1X < J < -. (6.7) 

^ 2 = H^ 2 + i) 5 " 2 

which is also the intersection of A = with S = in this case. Hence there is no stationary 
point in the open region Q 33 when 1 < 7 < 2; all the stationary points are on the boundary 
of 33 in this case. 




Fig. 6. 1 . The domains Q. 3 and Q. ^ 

We claim that solutions inside ^33 do not leave ^33 from its sides (excluding possibly 
edges, corners, or stationary points) as s increases for all 7 > 1. Note that the sides of ^33 
in the surfaces K = 0, or J = or H = are invariant regions. We need only to prove 
that no solution leaves O33 from the top side K = 1 — I. The top side has an outward 
normal 

n= (1,0,1) 
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arranged in the coordinate order (I, J, K) . We calculate the inner product of the normal 
n with the tangent vector of an integral curve of (6.1-3) on the top surface to yield 

n--^(I, J, K) = (l- I)(B + \KA) < 

which proves the claim, since > in fi 33 . 

We analyze the local structure of solutions at the stationary point (6.6). Let 

7=/ _ 7 , J = J, K = K-1^. 
To first order system (6.1-3) can be reduced to 

dl - ^-^a + V2K) 



dr 7 3 

dl _ (2- 7 )( 7 -l) 2 
dr 2 7 3 



J 



dK 



dr v/2 7 3 v 2 

Its three eigenvalues are found to be 

(7 -I) 2 



A± = 



A 2 = 



2 7 3 

(2-7)(7"l) 2 



7 + 1 ± V / (7 + l) 2 +4 7 ( 7 -l)] 



2 7 3 

Corresponding eigenvectors are found to be 

v± = (4A 0, 7 - 3 ± v / (7 + l) 2 +4 7 ( 7 -l)) 

v 2 = (0, 1, 0). 

It follows that A+ < 0, A2 < 0, and A_ > for 1 < 7 < 2. Hence the stationary point (6.6) 
is hyperbolic. The local structure is depicted in Fig. 6.2. 

To determine monotonicity of solutions, we need to deal with the surfaces A = and 
B = inside 33 . Consider first the surface 

B = 1(1- I) 2 + (1- I)J 2 -2IK 2 = in < I < 1. (6.8) 
An outward normal is given by in the coordinate order (/, J, K) 

(2K 2 - (1 - I) 2 + 21(1 -I) + J 2 , -2(1 - I) J, 4IK) =n B . 
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We calculate the inner product of the normal n b with the tangent vector of an integral 
curve of (6.1-3) on the surface (6.8) to yield 

n B -^(I,J,K) = (l-I) 2 HJ b , 

where 

J b = ( 7 /-l)(l-/)-(2- 7 )J 2 . 

See figure 6.1 for the position of J b = 0. So integral curves of (6.1-3) go from B > to 
B < on the surface B = with J b < 0, but reverse their directions on the surface B = 
with J b > as r increases. 

I K 




1/2 1/y 

Fig. 6.2. Local structure of stationary point (6.6). 
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Now consider the surface 

A = 2(1- I) 2 + (7 - 1)[J 2 - 1(1 - I)} - 2K 2 = 0, < I < 1. 
An outward normal is given by 

(4(1 - I) + (7 - 1)(1 - 2/), -2(7 - 1) J, 4iT) = n A 

in the order (J, J, K). We similarly calculate the inner product of the normal ha with the 
tangent of the integral curves on this surface to yield 

n A ~(I,J,K) = HJ a , 

where 

J a = [7 + 3 - 2(7 + 1)/](1 - 7)(1 - 7O + (7 " 2 </ 2 (l " 20- 

See figure 6.1 for the position of J a = 0. In particular the cylindrical surface J a = 
intersects the (I, 0, iO-plane at I = I/7 and the surfce H = at I = 2 (lj-i) • Therefore 
integral curves of (6.1-3) go from A > to A < on the surface A = with J a < 0, 
but reverse their directions on the surface with J a > as r increases. The portion of the 
surface A = below K = will not be used. 




We need to study the local structure of solutions near stationary points given in 
formula (6.5). It is helpful to first describe solutions of (6.1-4) with data u$ = 0, vq > 
0, po > 0. Our far-field solutions (3.4-5) can be written in terms of 7, J, K as follows: 



I = s 2 v 2 




s 2 v 2 (6.9) 
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valid for s G [0, s*] where s* = 1/r*. These solutions are all in the surface H = and they 
all start from the origin (J, J, K) = (0, 0, 0) and end at points of the stationary edge given 
by (6.5). As the initial Mach number Mq = |i>o|/ c o varies in (0, oo), the ending points of 
the solutions cover all the interior points of the stationary edge (6.5) exactly once. 

Now we study the structure of solutions near the stationary edge (6.5). More precisely, 
let us use E to denote the open set 



E = {{I, J, K) = (a, y/a(l -a), I - a) | Va € (0, 1)} (6.10) 

which contains the interior points of the curve (6.5). At any point of E, the linear part of 
the right-hand side of (6.1-3) is found to be given by the matrix 



-(1-a) 



(l-a)(l + 2a) -2(1- a) J a 4a(l - a) 

-2(2a-l)J a -2(2a-l)J, 

(1-a) p=±(l -2a) + 2(1 -a)] -(1 - a) (7 - 1) J a 2(1 - a) 



ct 

2 



where J a = a/ a(l — a). The eigenvalues and eigenvectors of this matrix can be found 
as follows. Since the whole line E is stationary, we expect that Ai = is an eigenvalue, 
and the tangent vectors of E are associated eigenvectors. The solutions of (6.9) offer 
another set of eigenvectors which can further simplify the calculation. We find the other 
two eigenvalues to be 

A 2 = -2(1 - a) 2 (l + a) < 0, A 3 = -(1 - a) 2 (l - 2a) 

with associated eigenvectors 

1 — 2a 1 — a 

f2 = (2Va, 7= , —1=-) 

VI — a v a 

n 3 = (2(l-a)[2(7-l)a-l], -2( 7 + 1)^(1 - a)(2a - 1), (6.11) 

-(l-a)[2(37-l)a-(7 + 3)]). 
Along this direction (6.11), we calculate the following information: 

n 3 • n A \ E = 2(l-«)(2«-l)[2(3 7 -l)«-( 7 + 3)] (6.12) 
n 3 -n B \ E = 2(1 - af{\ - 2a)[2( 7 - l)a - 1] (6.13) 
n 3 -n H \ E = 2(1 - a)(l - 2a)(l + 4a) (6.14) 

where is used to denote an outward normal to the surface J 2 = 1(1 — I). 

From now on we restrict ourselves to the case 1 < 7 < 2, since we have done the 
other cases 7 > 2 to earlier papers. The integral curves along n 3 starting from E in the 
direction of increasing s are depicted in Fig. 6.3. 
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It can be seen that n 3 points into different directions relative to the positions of A = 0, 
B = 0, and H = and the sign of the eigenvalue A3 changes when a varies. Our goal is to 
use the appropriate direction to construct solutions with monotone increasing s. Relative 

to the surface B = we find that n 3 points into fi 3 for a E ( 2(7-1) ' ^)' ^ n ^° ^ 33 an< ^ 

B < for a E (|, 2(7-1) ) ' anc ^ — 77-3 P°i nts i nto ^3 f°r ^ e (0> §)• Relative to the surface 

A = we find that n 3 points into 33 and A < for a E ( 2(37^-1) ' -*-)' ^ nto ^. > for 

a G (5, 2(37"- 1) ) ' anc ^ — n 3 points into O3 for a E (0, |). 

Now we can depict the integral curves inside Q33, see Figure 6.4. 




(Yf3)/(&y-2) 



Fig. 6.4. Integral curves in the intermediate field. 
First we observe that s is an increasing function of r inside n 33 . J is an increasing 
function of r if < I < |, but changes to decreasing when I E (|, l). It is useful to 
observe the stable manifold of the system (6.1-3) at the point (6.6) which contains the 
transitional integral curve in the case vq = considered in Sect. 5 and the transitional 
integral curve from a point of E with a critical I = Ih- There are three kinds of integral 
curves depending on their positions relative to the stable manifold. The first kind consists 
of integral curves which are below the manifold and go to the stationary point (1, 0, 0). 
Each of the second kind is right on the manifold and goes to the stationary point (6.6). 
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Each of the third kind is above the manifold and goes to a stationary point on the curve 
E with 7 G (0, |) given in (6.10). We mention in particular that no integral curve from 
inside O33 goes to a point of E between | < 7 < 1 because 7 is an increasing function of 

r and 113 is pointing towards (1, 0, 0) for < I < 1. See Lemma A.l of the Appendix for 
a complete proof. Also there is no integral curve from inside ^33 which goes to the point 
(0, 0, 1) because J is an increasing function of r for 7 G (0, |). 

We need to extend the intermediate field solutions to all s. Integral curves ending 
at point (6.6) are already defined for all s > since r — > +00 and the right-hand side of 
equation (6.4) does not vanish at (6.6). For the first kind of integral curves which end at 
(1, 0, 0), we show that s approaches finite values although r — > 00. Since the proof is long 
and tedious, we put it in the Appendix to avoid interuption of our construction. We use 
the natural continuation of vacuum p = to extend our solutions till s = 00. The values 
of u, v are not needed to be specified in the vacuum. Near the edge of the vacuum we find 

**\/3 a — 7 + 1 



v = b{r-r**) p , (3 



2( 7 -l) 



where b > is any constant. The asymptotics of u and c given in Sect. 5 are still valid. 
From these asymptotics we find that M = 00 and M s = , which are the same as in 
Sect. 5. 

For the third kind of integral curves which end on the upper half of E with a G (0, |), 
we show that s approaches finite values also. The proof is easy because we have 



ds s 1 — 7 



dJ J 1-27 



(6.15) 



from the second equation of system (4.3). The right-hand side of (6.15) is nonsingular for 
(I, J,K) G E with < I < \. Thus s is finite around any point of E with < I < \. 

We comment that integral curves on the surface J 2 = 7(1 — J) starting from the 
origin and ending on a point of E with a G (7^, 1), where 7^ G (I/7, 1/(2(7 — 1))? can b e 
continued through the direction (6.11) into 33 , and they go toward (1,0,0) with finite 
ending s values. The further extension by vacuum is also valid. But integral curves on 
the surface J 2 = 7(1 — 7) starting from the origin and ending on a point of E with 
a G (1/2, 7^) are continued along the direction (6.11) into ^33 and curl back to a point 
of E with < a < 1/2. Finally integral curves on the surface J 2 = 7(1 — 7) starting 
from the origin and ending on a point of E with < a < 1/2 will be continued along 
direction (6.11) into 77 > and discussed in the next section. It would be interesting to 
find explicitly the function 7^ = Ih(l)- We know that 7^ = 1/2 for 7 = 2, see [16, 17]. 
Numerically, we find that 

Zfc(l.l) = 0.93, I h (lA) = 0.77, I h (1.7) = 0.63, 

see paper [16] for numerical procedures used. For the special data (uo, vq, po) with w = 0, 
the critical value 7^ can be expressed through the solution formula (6.9) by the critical 
initial Mach number 

M h ( 7 ) = y/h/(l-I h ). (6.16) 
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Since I h G (I/7, 1/(2(7 - 1))), we find that 



^ < Mh (j) < (6.17) 



7-1 X " 27-3 

The lower bound is valid for all 7 G (1, 2), but the upper bound is valid only for 7 G (3/2, 2). 
So we have constructed all solutions globally except those which end on E with < I < \. 

7. Inner-field solutions. 

We extend solutions that end on the set E of (6.10) with < / < | in this section. 

We find that these solutions go along the directions n% given in (6.11) into the region 
between A = and S = 0in0</<^ and J 2 > 1(1 — I), and eventually go to infinity. 
The scaled variables /, J, and K are not suitable for this portion of the solutions. 

We restart from system (2.7) with data (u, v, p, r) satisfying the relations E : K = 1 — 1 
and J 2 = 1(1 — I) for / G (0, |). In terms of (it, v, c, r), these data are in the form 

u\ =13 

I r=a ' 

v \r=a = V/3(« - J) (7-1) 

c| = a — 3 

I r=a 



where a > and (3 G (0, ^) are arbitrary. The directions 77-3 given in (6.11) point into the 
region A > 0, > 0, E > 0. 

Asymptotic analysis shows that problem (2.7) and (7.1) has solutions (-u(r), v(r),p(r)) 
which vanish as r — > 0+. The solutions are in the region A > 0, > 0, E > 0. We perform 
an asymptotic analysis a priori to determine the orders at which (u, v, p) vanish as r — > 0+. 
For polytropic gases p(p) = A 2 /9 7 , we find 

ar k d 



where 



2(2-7) 2-7 



and a > and d > are arbitrary. We can easily verify that these directions are such that 
A > 0, 6 > 0, E > for r > 0. In terms of the sound speed c, we find 

= fcy/A 2 (2-7) 

° (lnf)(7-i)« ' 
We are motivated by this asymptotic analysis to use the scaled variables 

U=-,V=-,R=-. (7.3) 

c c c 
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We can rewrite system (2.7) into a new form 



f du r— ii 

dr ~ c 3 

dv — uv /\ 
dr "c 3 " 



dr 
dr 



dc 7— 1 r—u r\ 

Wt" 2 



(7.4) 



where r is a parameter. In terms of the variables in (7.3), we find 

|Z = (R - u)t - XU(R - U)Q = A t 
% = UVA - XV (R - U)Q = d 
^§ = R(R - U)A - XR(R - U)Q = B x 



(7.5) 



where A = and A = 1 - (U - R) 2 , 6 = V 2 - U(R - U) and S = (R - U)0 - UA. 
System (7.5) is autonomous for (U, V,R). We find that the last equation in (7.4) can be 
written as 

dc 

-- = X(R - t/)e. (7.6) 
car 

So c can be integrated from (7.6) once (U, V, R) are obtained from (7.5). The corresponding 
data of (7.1) for (7.5) and (7.6) are any stationary point (U*, V*,R*,c*) satisfying 



R* -U* = 1, V* =U\ < 2U* < R\ 



c* > 0. 



The asymptotic analysis (7.2) shows that 

(U, V, R) -> (0, 0, 0), r^0 + 



(7.7) 
(7.8) 



(7.9) 



After (7.5-9) are solved, we use the third equation in (7.4) to show that r is an increasing 
function of r G R and r — > or a G (0, +00) as r goes to =Foo respectively. 

For more refined change of variables, we further find that our asymptotic analysis 
(7.2) implies that 



(7\(7-l)«-l 

U = abr In — 



(7.10) 



R = br 



fin- 



( 7 -i)a 

r / 
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where we set b = (ky/(2 — r ))A 2 ) 1 . It follows that the above asymptotic behavior yields 
the function relations of {R, U) to V 

( b 2 k 2 \ fbk\ 2(l - 1)a 

R= 7^ fexp v^J \v) 

U = -RV 2 
2 

which are not analytic at V = 0. It can be shown that the directions (7.10) enter the 
region A\ > 0, B\ > 0, and C\ > 0. We omit the proof. 

We transform the direction 773 of (6.11) to the new variables (U, V, R) at a point of 
(7.7). It is possible to use linearization for the derivation, but we prefer a direct transfor- 
mation. Let n= (ni, 77-2, 77.3) be an arbitrary vector in the coordinate space (I, J, K), 
which is a tangent direction of an integral curve of (6.1-3). So we have 

dJ 772 dK 773 

dl n\ ' dl n\ 

along the integral curve. Because the variables (I, J, K) are related to (U, V, R) by 

/ = U/R, J = V/R, K = 1/R, 
we find that the differentials (dR, dU, dV) satisfy 

RdV - VdR _ 77 2 -dR _ n 3 

RdU-UdR ~ tT? RdU -UdR ~ th' 

which can be solved to yield 

dU Uns — ?7 1 e/V ^773 — 772 
(ii? i??73 ' dR -R77.3 

Thus the corresponding tangent vector in the new coordinate space (U, V, R) of n = 
(ni, ?72, 773) in the coordinate space (I, J, K) is a multiple of 

(E/713-rai, F?73-?72, i?n 3 ), (7.11) 

where (ni, 772, 773) as well as (t/, V, i?) are evaluated at the same point (u, v, c, r). 
For our particular problem 773 of (6.11), we use a = U/R for a point of (7.7) to find the 
corresponding direction to be a multiple of 

n* 3 =(2 + 3(3 - >y)U - 5( 7 - l)U\ VU[(7 - &y)U - (7 - 1)], 

(l + t/)[ 7 + 3-5( 7 -l)C/]). (7.12) 
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We determine the relative positions of the regions where Ai > 0, B\ > 0, or C\ > 
in the region R > U, U > 0, V > 0. The common intersection A\ = Bi = C\ = is the 
curve (7.7) of stationary points for all < U < R. The intersection A\ = Bi = is in the 
plane R = jU. The intersection B\ = C\ = is in the plane R = 2U. The intersection 
A\ = C\ = is in the plane R = jz^U ' . We consider two domains 

On ={(U, V, R)\U>0, V>0, R> 2U, A 1 > 0, B 1 > 0} 
42 = {(U, V, R) | U > 0, V > 0, 2C/ > R > jU, A 1 > 0, d > 0}. 
See Fig. 7.1. 



V 




Fig. 7. 1. £2 41 and £2 
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We show that integral curves of (7.5) on the surface A\ = of the boundaries of both 
4 i and 42 are inward. In fact, the surface A\ = is represented by 



V 



2 U[2-(-y-l)U(R-U)] 



with an outward normal 



n A \= n u -2V, 



2R- (rf + l)U 

U[4 - (7 - 1) 2 U 2 } 



[2R- (-y+l)U]< 



where n\ is some expression that we do not need to compute. We calculate the inner 
product of this normal with the tangent direction of integral curves of (7.5) 

±iU, V, R) . n A1 = -2VC, - ^±^m Bl < 

since both C\ > and B t > on Ai = 0, in the range R > . 

We show next that integral curves of (7.5) are incoming on the surface Bi = of the 
domain f^i where R > 2U. In fact, the surface Si = is represented by 

XV 2 = 1 + XU(R -U)-(R- U) 2 

with an inward normal 

n m = (X(R-2U) + 2(R-U), -2XV, XU - 2{R - U)) . 
We find the inner product 
d 



dr 



(U, V, R)- n m = Ai[X(R - 2U) + 2(R - U)] - 2XVd > 



since Ai > and C\ < in the region R > 2U on the surface B t = 0. 

We show next that integral curves of (7.5) are incoming on the surface C\ = of the 
domain where 2U > R > ■yU. In fact, the surface C\ = is represented by 



XV 



U 



R-U 



- (1 - X)U(R-U) 



with an inward normal 

R 

ncv 



(R-uy 



+ (1 - X)(2U-R), -2XV, -(1-X)U- 



U 



(R-uy 



We find 
^(U, V, R)-n C i = 



R 



(R-U) 2 



+ {1-X){2U-R) 



A 1 



{l-\)U + 



u 



(R-U) 2 



B, 



> 
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since A\ > 0, B t < on C\ = in the region g^-t^ < R < 2U which contains -fU < R < 
2U. 

It is easy to see that integral curves of (7.5) are incoming on the surface U = of the 
domain Q 4 i since 

^ = R 2 V 2 > 
dr 

on U = 0. 

We show that integral curves of (7.5) go from fi 41 to fi 42 on the common interface 
R = 2U . A normal of R = 2U pointing from fi 41 to fi 42 is n= (2, 0, —1) and 

4-(U, V, R)- n= 2U 2 (1 -U 2 - V 2 ). 
dr 

on R = 2U. Using Ai > in Q 4 i and 42 on R = 2U, we find 

V 2 + U 2 - 1 > — > 0, 

1 — A 

thus integral curves of (7.5) go from 4 i to 42 on the interface R = 2U . 

So integral curves in 4 i exit 4 i only through the side R = 2U or the stationary 
edge (7.7). 

We need to study the local structure of solutions of (7.5) at the origin (U, V, R) = 
(0, 0, 0). In particular, we would like to prove the asymptotic behavior (7.2). To reduce 
the order of degeneracy, we introduce the variables 

X=-, V 2 = V 2 , R 2 = R 2 . 

JrC 

We find that 

X = a(ln^ , V 2 =2X (7.13) 

as r — > 0+ from the asymptotic analysis. So the origin (X, V 2 , R 2 ) = (0, 0, 0) will still 
be the place for us to look for outgoing solutions. Introducing t' by 

dr = Rdr, (7.14) 

we find the equations (7.5) in terms of (X, V 2 , R 2 ) to be 



ax 

dr 



: 7 = (1-X) [V 2 (l -X)-2X + X(l - X) 2 R 2 ] = (1 - X)A U (7.15) 



^ = 2V 2 {X - X(l - X) 2 R 2 - A(l - X) [V 2 - X(l - X)R 2 }} (7.16) 
^ = 2R 2 (1 - X) {1 - (1 - X) 2 R 2 -\{V 2 - X(l - X)R 2 ]} (7.17) 
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where we have set the expression in the brackets in (7.15) to be An. We need the local 
structure of (7.15-17) at (0,0,0). Linearization at the point (X, V 2 , R 2 ) = (0,0,0) yields 



f^ = -2X + V 2 



dV 2 _ Q 



dr 



dR 2 
dr' 



which has eigenvalues 



2R 2 

Ai = —2, A2 = 2, A3 = 0, 



and associated eigenvectors 

(X, V 2 , R2) = (1,0,0), (0,0,1), (1,2,0). 




Points of (7.7) 



Fig. 7. 2. The center-unstable manifold of the origin. 

By the Center Manifold Theorem [8], system (7.15-17) has a C k center manifold which 
is tangent to the vector (1,2,0) for any k < 00. Suppose the center manifold takes the 
form 



X = g(V 2 ), 0(0) =0, </(0) = l/2 
R2 = h(V 2 ), h(0) = h'(0) = 0. 



(7.18) 
(7.19) 
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Then the flow on the center manifold is described by 



dV 2 



2V 2 {g - g(l - gfh - A(l - g) [V 2 - g(l - g)h\) . 



(7.20) 



dr' 



Using (7.18-19), we can approximate (7.20) to find 



dV 2 



(2- 7 )U 2 2 + 0(U 2 3 ). 



dr' 



So the center manifold is unstable since 2 — 7 > 0. Now we consider the center-unstable 
manifold consisting of the unstable manifold and all center manifolds, see Kelley's paper 
[8] for the existence and smoothness. This manifold is two-dimensional. It can be shown 
easily that they enter the region A\ > 0, B\ > 0, C\ > 0. Furthermore, we can use 
Henry's approximation of center manifold ([7,8,2]) to find 



which can be used to show that all center manifolds enter the region An > 0. 

The center- unstable manifold enter the region 41 . Note that the unstable manifold, 
i.e., the R 2 axis, ends on one end of (7.7). We conclude that nearby integral curves in the 
center-unstable manifold end on (7.7) too, see Fig. 7.2. 

By our earlier analysis, we know that integral curves in the center-unstable manifold 
can only exit O41 through the divider between O41 and O42 or points on (7.7). From the 
continuity of the center- unstable manifold, we conclude that every point on (7.7) is covered 
by the center-unstable manifold. 

Finally, we show that r — > 0+ at the origin (X, V 2 , R 2 ) = (0, 0, 0). It is easy to see that 
c is an increasing function of r from (7.6) since O is positive along the center manifolds. 
Thus c is bounded as r — > —00. From (7.3), we have r = Rc. Hence r — > since R — > as 

T — > —OO. 

We point o ut that both th e Mach number M = \/U 2 + V 2 and the pseudo Mach 
number M s = y/ (U - R) 2 + V 2 go to as r — > 0. Also, both U and R — U go to faster 
than V, which explains the spiralling phenomena in both the physical space and in the 
self-similar coordinates. 

Thus each point on (7.7) has a solution going to the origin through O41 with r — > 0+. 

8. Conclusions. 

We summarize our results. By a weak solution to the 2 — D Euler equations (2.1) 
we mean a bounded vector function (u, v, p) satisfying the equations in the sense of dis- 
tributions. Since we deal with only continuous weak solutions in this paper, we shall not 
mention any other requirements such as Rankine-Hugoniot relation or entropy conditions, 
let c = \/p'(p) be the sound speed. Let M = \/ u 2 + v 2 /c be the Mach number. Let 
M s = [(u — £) 2 + (v — r]) 2 ] 1 / 2 jc be the pseudo Mach number of the self-similar flow in the 



9=lv 2 -^V 2 2 +0(Vi) 



h 



o(v 2 % 
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self-similar coordinates (£,rj). In the notation of system (2.7), the pseudo Mach number 
M s = [{u-rf + v 2 ] 1 / 2 /c. Set 



M o = y «o + v o/ c o, c o = Vp'(Po) 

which are in consistency with previous versions of M (Sect. 5 & 6). Also, let us introduce 
the vector parameter ( "Mach vector" ) 

V c o c J 

Recall that we can assume vq > without loss of generality. Also, the case vq = has 
been resolved in [17]. 



| u 

V 

p 

p 

n 




r o 




— v o 

V 

u 

- 



Or x 



Fig. 8.1. 



A solution {u, v,p) vs. the x-axis at time. 
t = 1 with a datum uq = and M < v2. 

Theorem. Assume p(p) = A^p 1 for any constants A2 > and 1 < 7 < 2. Then for any 
datum (uq,vq,po) with uo > 0, vo > 0, and po > 0, there exists a weak solution (u,v,p) 
to the initial value problem (2.1) and (2.5). The solution is continuous, self-similar, and 
axisymmetric for t > 0. It takes on initial datum almost everywhere and in L^ oc for any 
q > 1 as t — > 0+. Additionally, we have the following properties. 

(I). If uo = 0, the solution consists of two or three smooth pieces, depending on the 
initial Mach number Mq. The first piece of the solution is always given explicitly by the 
formula (3.4) with constant Mach number near r = 00. 

(La). If M < \pl, the solution consists of two smooth pieces. The end point of 
the first piece of the solution is on the set of stationary points (6.5) with an I < 1/2, 
which is transformed to be the stationary set (7. 7) from which the second piece begins in 
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the direction given in (6.11) or equivalently (7.12). It is an integral curve on the center- 
unstable manifold of the ODE (7.15-17) at the origin. See Fig. 7.2 and 8.1. Both the 
Mach and the pseudo Mach numbers are zero at the center point of the solution. 

(Lb). If M > \f2, but M < M/j( 7 ) (which is defined in (6.16) ), the solution consists 
of three smooth pieces. The second piece is given by the solution of the ODE(4.3) with 
direction (6.11) such that the solution goes to a stationary point on (6.5) with an i" < 1/2. 
The third piece is identical to the second piece of a solution described in the previous case 
Mo < V2. 

(I.c). If Mo = M/j( 7 ), the solution consists of two pieces and the second piece is given 
by the solution of the ODE(4.3) with direction (6.11) such that the solution goes to the 
stationary point (6.6). 

(I.d). If Mo > M/ l (7), the solution consists of three smooth pieces and the second 
piece is given by the ODE (4.3) with direction (6.11), see Fig. 8.2. The radial and pure 
rotational conponents and the density function of the solution are increasing functions of 
the spatial radius r in the region u**t < r < r*t where r* is given by (3.5) and u** > 
is the radial velocity at the inner end of the second piece. The third piece is the vacuum 
p = with domain < r < u**t. At the edge of the vacuum, there is no rotation, the 
Mach number is infinity, and the pseudo Mach number M s = (for 7 < 3). 

(I.e). The critical Mach number M/ l ( 7 ) has a lower bound 

V7/( 7 -l)<M h ( 7 ) 
for all 7 G (1, 2). It has an upper bound 

M h ( 7 )< V2( 7 -l)/(2 7 -3) 

if 7 G (3/2,2). 

(II). Assume uq > 0, po > 0, and vq > 0. Depending on the vector parameter Mo, 
the solution can be globally smooth, contain a region of vacuum as in case (La), or contain 
an inner piece of a solution of (I.d). In terms of the selfsimilar coordinates and variables, 
the ODE (6.1-3) has at the point (6.6) a stable manifold, which contains the transitional 

solutions in the special cases uq = of case (I) and vo = 0. If the vector parameter M 
is such that the solution is on the stable manifold, then the entire solution consists of one 

smooth piece and is governed by the system of ODE (3.1). If the vector parameter M is 
such that the solution is below the stable manifold, then the solution consists of a vacuum 
region at the center and the far-field solution given by the system of ODE (3.1) with data 

(3.2). If the vector parameter Mq is such that the solution is above the stable manifold, 
then the solution consists of two pieces, one is given by the system of ODE (3.1) with data 
(3.2), the other by an integral curve in the center-unstable manifold near the origin of 
(7.15-17), and the density p of the solution vanishes only at one point, the origin of space. 

We omit the proof of the theorem with the remark that it can be verified that the 
pieces of integral curves when combined continuously in the order specified in the theorem 
form weak solutions to (2.1)(2.4). 
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The case 7 = 1. We now point out necessary changes in Sections 4-7 to accommodate 
the case 7=1. The equations (4.3) are still valid although the equation for K is trivial 
since the sound speed c is a constant. We need to supplement (4.3) with the equation for 
p in (3.1), the equation for p being decoupled from (4.3) now. Conclusions in Sect. 5 are 
modified as follows: all integral curves in Q go to the point (I, K) = (0, 1) except the trivial 
one K = 0. In Sect. 6 all integral curves in 3 = fi 33 from the origin (/, J, K) = (0, 0, 0) 
enter the set E with a < 1/2. In Sect. 7 the conclusion is the same, but the proof, which 
we shall omit, needs nontrivial changes since the domains O41 and are not bounded 
anymore. The overall conclusion for 7 = 1 is the same as stated in the previous theorem, 
the main difference is that the threshold Mach number Mh(l) = 00 and the threshold 
stable manifold mentioned in Step II is the bottom surface K = 0, thereby all solutions in 
the case 7 = 1 have point cavities at the centers. 




Fig. 8.2. 



A solution (u, v, p) vs. the x-axis at time 
t = 1 with a datum uq = and M > Mh("f). 
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Appendix. Finiteness of the parameter s at point (1, 0, 0) 

In this appendix we show mainly that the parameter s, the reciprocal of the radius r, 
approaches a finite value as an integral curve approaches the stationary point (1,0,0) in 
^3 defined in section 6. 

Lemma A.l. Assume 7 > 1. Then for any integral curve of the system (6.1-4) that goes 
from I < 1/2 to I > 1/2 inside the domain O33, there exists an e E (0, 1) such that the 
integral curve is inside the cylinder 

C e = J 2 -el(l -I) < (Al) 

whenever its I component is in (|, l). 

Proof. We choose e G (0, 1) such that the integral curve is inside the cylinder J 2 < 
el(l — I) at the special point I = ^. For this e, we show that integral curves on the 
surface of the cylinder J 2 = el(l — I) inside O33 in the portion I G (|, l) are all going 
into the cylinder J 2 < el(l — I). We first calculate an outward normal to the cylinder 
J 2 = el(l - I) to be 

n £ = (-e(l-2/),2J,0). 

We then calculate the inner product of the normal n £ with the tangent direction of any 
integral curve on the cylinder J 2 = sl(l - I) 

n £ - A (j, j, K)\ J2=£I{1 _ I} = -e(l - 2J)(1 - J)[J(1 - I) 2 (l + e) - 2IK 2 ] 

+ 2el(l - J)(l - 2/) • [(1 - J) 2 - K 2 ] 
= £/(l - 27)(1 - + e)(l - I) 2 + 2K 2 + 2(1 - I) 2 - 2K 2 } 

= e(l - e)I(l - I) 3 (l - 21) < for I e Q, 1^ . 

Hence the chosen integral curve will remain inside the cylinder J 2 < el(l — l) in / e l). 
This completes the proof of Lemma A.l. 

Remark. From the proof of Lemma A.l it is clear that any integral curve entering 33 
from the edge E (6.10) will get into a cylinder J 2 < el(l — I) for some e G (0, 1) for / 
close to 1. 

Lemma A. 2. Assume 7 > 3/2. Then for any integral curve that ends at (1, 0, 0) of the 
system (6.1-4) in the domain Q 3 , there are three numbers e G (0,1), j3 G (2, 00), and 
I G (|, l) such that the integral curve is inside the cylinder C e < and beiow the surface 

Bp = (1- I)J 2 + 1(1- I) 2 - /3IK 2 =0 (A2) 

for all I e (1,1). 
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Proof. We first compute an outward normal rip to the surface Bp = 0. 

np= ((3K 2 - (1 - I) 2 + 21(1 + J 2 , -2(1 - I) J, 2(3IK). 

We calculate the inner product of this np with the tangent vector of any integral curve on 
the surface Bp = 0: 

= [f3K 2 - (1 - I) 2 + 21(1 -I) + J 2 ] • (1 - i") [(1 - I) 2 1 + (1 - I) 2 - 2IK 2 ] 
- 2(1 - I)J 2 (1 - 2J)[(1 - I) 2 - K 2 ] + 

PIK 2 (1 - I) [( 7 - 1) J 2 - ( 7 - 1)7(1 - I) + 2(1 - I) 2 - 2K 2 ] 
21(1 - /) + i J 2 ] (1 - I) 2 (l - |) [J 2 + 1(1 - I)] 

+ 2(1 - I) 2 (l - 2I)J 2 ^ [J 2 -((3- 1)1(1 - I)] + (1 - I) 2 [J 2 + 1(1 - I)} j- 



(7-1)/ -_(!-/) 



J 2 -(7-lK 2 (l-/) + 2(l-l)/(l-/) 2 



Rewriting the last expression in the form of a polynomial of J, we obtain 



^{>-/) 2 (l-|) + ^(l-/) 2 (l-2/) + ^ 

+ J 2 {2/(1 - If (l - |) + (1 - ') 3 (l - |) - 2^(1 - /) 3 (1 - 21) 

+ 2 ( 1 -^) (i-^) 4 -(7-iK(i-/) 3 } 



( 7 -!)/--(!-/) 



( 7 -!)/--(!-/) 



+ 2/ 2 (l - if ( 1 - |) - (7 - 1)/ 2 (1 - If + 2 (l - I ) /(I - /) 5 



= J 4 (1 . /)2 {/3 - 2 + [0( 7 - 1) - 2]/} + +J 2 (1 - /) 3 



PI 

+ 1(1 -If 



2/(l-|)-(7-l)/ + 2 (l-I)(l-7) 
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Observing that both coefficients of J 4 and J 2 are positive when f3 > 2, 7 > 2, and 
I E (|, l), we find in this case in the region C e < that 



+ el(l-l) 4 
+ 1(1 -I) 4 



4/U 4) +1 -| 



.--(,-1)7 + 2/!-- 



(-4-3) 



l-| + ( 7 -l-||I 



+ 4eJ 1 - 







+ 



+ e -£-(f + 7-l)/ + 2ll 



= /(l-/) 4 F( £ ,ft/, 7 ) 



5)} 



where F denotes the expression in the curly brackets. 

We show that the inequality (A. 3) still holds even when 7 G [3/2, 2). The quadratic 
polynomial of J 2 attains its maximum at 



r2 . _ 7(l-/)[4/(/?-l) + /?-4] 



At /3 = 2, we calculate 



{J )max 



2/[2-/3( 7 -l)] +4-2/3 



1-7 

2(2-7) 



+ 0((1-/) 2 ) 



(A4) 



(A5) 



is 



for / close to 1. Since (J 2 ) max is a continuous function of /3 near 2, we see that ( J 2 ) 7 
positive for all 7 G (1, 2) when (3 is close to 2 and I < 1 is close to 1. Hence the quadratic 
polynomial of J 2 is an increasing function of J 2 in the interval [0, (J 2 ) raa J. We show next 
that 

el(l -I)< (J 2 ) max (A6) 

for all 7 G [3/2, 2), and e < 1 when / is close to 1 and (3 is close to 2. In fact, (A. 6) holds 
by continuity, the fact (A. 5) and the inequality 



e < 



1 



2(2-7) 



which holds for all 7 G [3/2, 2). So inequality (A. 3) still holds in the region C e < even 
when 7 G [3/2,2). 
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We need F to be negative near 7 = 1 for some e > 0, j3 > 2, and 7 > 3/2. For any 
fixed integral curve that goes to the point (1,0,0), we first choose e G (0,1) so that the 
integral curve lies inside C e < for 7 G (7, 1) for some I < 1. Fix this e. We calculate the 
value of F at the extreme point (3 = 2, 7 = 1: 

F( £ ,2,l,7)= £ 2 ( 7 -2)+£- 7 + l 

= (l- e )[l + 2e-7(l + e)] 

<(l- e )[l + 2e-|(l + e)] 

= -^(i-^) 2 <o. 

Since F is a continuous function of (/?, J) near the point (2, 1), we conclude that for each 
7 > 3/2 there exists an I G (\, l) such that F(e,(3, 1,7) < for all 7 G [7, 1) when (3 is 
close to 2. Since 7?g = is close to B = 0, we can choose a /3 > 2 such that the integral 
curve lies below Bp — at 7 = 7. This integral curve will remain under the surface Bp = 
for all 7 > 7 because of the sign F < 0. This completes the proof of Lemma A. 2. 

Lemma A. 2 does not hold for 7 < 3/2. We need to deal with the case 7 < 3/2 
separately. Observe that the intersection point of the curve = with the stationary 
edge (6.5) is topologically different for 7 < 3/2, see Fig. 6.1. 

Lemma A. 3. Assume 1 < 7 < 3/2. Then integral curves in fi 3 of system (6.1-3) on the 
surface of the domain 

(2- 7 )J 2 < £ ( 7 7-l)(l-7), (A7) 

where e G (0, 1), always enter the domain. 
Proof. We calculate an outward normal to (A. 7): 

n= (-5(7 + 1-277), 2(2-7)7, 0). 

We then calculate the inner product of this normal with the tangent vector of an integral 
curve of (6.1-3) on the surface of the domain (A. 7): 

n 7, K) = e(l - 7)(2 7 7 - 7 - 1)5 - 2(2 - 7 )7 2 (27 - 1) [(1 - 7) 2 - K 2 ] (A.8) 

= e(l - 7) 2 (2 7 7 - 7 - 1) [7 2 + 7(1 - 7)] - 2(2 - 7 )7 2 (27 - 1)(1 - 7) 2 
+ K 2 [2(2 - 7)</ 2 (27 - 1) - 2e7(l - 7)(2 7 7 - 7 - 1)] 

= e(l - 7) 2 (2 7 7 - 7 - 1)[7 2 + 7(1 - 7)] - 2(2 - 7 )7 2 (27 - 1)(1 - 7) 2 
+ 2e(l ~I) 2 K 2 . 
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We observe that this inner product is an increasing function of K 2 . In the domain Q 3 , the 
largest K 2 is achieved on the surface B — 0. From the step (A. 8), we find that this inner 
product in negative for I > 1/2. This completes the proof of Lemma A. 3. 

Lemma A. 4. Assume 1 < 7 < 3/2. Then for any integral curve that gets into (at least 
once) the domain 

( 2 _ 7 )j2<( 7 j _i)(i_/) (A9) 

and ends at (1,0,0) of the system (6.1-4) in the domain fi 3 , there are three numbers 
e G (0, 1), (3 G (2, 00), and I G (|, l) such that the integral curve is inside the cylinder 
(A. 7) and below the surface (A. 2). 



Bp = (1-I)J' A + I(l-If - (3IK 2 =0 (A2) 



when I G (J, 1). 



Proof. The proof parallels that of Lemma 2: use the surface of the domain (A. 7) in places 
where C £ = is used. 

Theorem A. For all 7 > 1, the parameter s is finite for any integral curve of (6.1-4) that 
goes to the point (1, 0, 0) inside O3. 

Proof, (i) Assume 7 > 3/2 or the integral curve gets into the domain (A. 7) at least once, 
hence the integral curve remains in Bp > for I close to 1. We find from the first equation 
of (4.3) that 

The right hand side of (A. 10) is a decreasing function of K 2 in fi 3 . So for I close to 1 we 
find 

f ((3-1)1(1- I) - J 2 f (13 - 1)1(1 - 1) - J 2 v ; 

The very last expression of (A. 11) is an increasing function of J 2 , so we further find by 
using J 2 > that 

ds ~ y f } (3-1 (3-1 

So it can only take a finite amount of s for I to reach 1. (ii) Ifl< 7 <3/2 and the integral 
curve in entirely outside the domain (A. 9) when / is close to 1, we have 

( 2 _ 7 )j2 > ( 7 / _i)(i_/) (A12) 
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(A13) 



for / close to 1. We can use equations (6.2) and (6.4) to find 

dJ_ _ J(f - 21) 
~ds ~ s(l-I) ' 

Using (A.12) in (A.13), we find 

dln(l) _2I-1 2 1 2 
dins 1-1 K J> 

where a 2 is a positive constant. It is easy to derive that s < oo from this last differential 
inequality. 
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Figure Captions. 



Fig. 1.1. The velocity vector of a typical swirling solution. Here 7 = 1.4, uo = 0.0, 
vq = 1-0, po = 1-0, and Ai = 1.0. The figure shows the solution at time t = 1 in the square 
|x| < 1.7535, and \y\ < 1.7535. 

Fig. 1.2. The fluid density of the typical swirling solution in Fig. 1.1. Outside the 
square the density is the constant 1. 

Fig. 5.1. Phase portrait of solutions without swirls. 

Fig. 5.2. Estimate of the transitional Mach number M/^7) < -^y. In this figure 
d = (7 - l) 2 /2, / = M" 2 ( 7 ), and N and D are where the numerator and denominator of 
the right-hand side of (5.8) vanish respectively. 

Fig. 5.3. A zero-swirl transitional solution. 

Fig. 5.4. The parameter s is finite near (1, 0). 

Fig. 6.1. The region Q3 and O33 . 

Fig. 6.2. Local structure of solutions at the stationary point (6.6). 
Fig. 6.3. Solutions in the intermediate field for data with u = 0. 
Fig. 6.4. Integral curves in the intermediate field. 
Fig. 7.1. The domains O41 and O42. 
Fig. 7.2. The center-unstable manifold. 

Fig. 8.1. A solution (u, v, p) vs. the x-axis at time t = 1 with a datum uq = and 
M < V2. 

Fig. 8.2. A solution (u, v, p) vs. the x-axis at time t = 1 with a datum uq = and 

M > M h {n). 
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Fig. 1.2 




